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ABSTRACT 

The presence of VI ions can be indicative of plasma temperatures of a few x 10 5 K that is expected 
in heat conduction layers between the hot shocked stellar wind gas at several 10 6 K and the cooler 
(10 4 K) nebular gas of planetary nebulae (PNe). We have used FUSE observations of PNe to search for 
nebular O VI emission or absorption as a diagnostic of conduction layer to ensure the presence of hot 
interior gas. Three PNe showing nebular VI, namely IC418, NGC2392, and NGC6826, have been 
selected for Chandra observations and diffuse X-ray emission is indeed detected in each of these PNe. 
Among the three, NGC 2392 has peculiarly high diffuse X-ray luminosity and plasma temperature 
compared with those expected from its stellar wind's mechanical luminosity and terminal velocity. 
The limited effects of heat conduction on the plasma temperature of a hot bubble at the low terminal 
velocity of the stellar wind of NGC 2392 may partially account for its high plasma temperature, but 
the high X-ray luminosity needs to be powered by processes other than the observed stellar wind, 
probably caused by the presence of an unseen binary companion of the CSPN of NGC 2392. We have 
compiled relevant information on the X-ray, stellar, and nebular properties of PNe with a bubble 
morphology and found that the expectations of bubble models including heat conduction compare 
favorably with the present X-ray observations of hot bubbles around H-rich CSPNe, but have notable 
discrepancies for those around H-poor [WR] CSPNe. We note that PNe with more massive central 
stars can produce hotter plasma and higher X-ray surface brightness inside central hot bubbles. 
Subject headings: planetary nebulae: general - planetary nebulae: individual (IC418, NGC 2392, 
NGC 6826) - stars: winds, outflows - X-rays: ISM 



1. INTRODUCTION 

Planetary nebulae (PNe) consist of the stellar mate- 
rial ejected by low- and intermediate-mass stars near 
the end of their evolution, before turning into white 
dwarfs. The physical structure of a PN is largely de- 
termined by the photoionization of the slow, dense wind 
ejected during the asymptotic giant branch (AGB) phase 
by the intense stellar radiation field and its interaction 
with the subsequent fast, tenuous wind emanating from 
the hot stellar core at terminal velocities up to ^4,000 



ture and raises the density of the hot gas ((Weaver et al.l 
Il977j ). significantly increasing the X-ray emissivity. Hy- 
drodynamic models of PNe with heat conduction predict 
diffuse X-ray emission that should be ea sily detectable 



with modern X-ray observatories Ce.g.. IVolk fc 



1985HZhekov fc PermottdTl998l:ISchonberner et al.l 
Steffen . Schon berner. fc Warmuth II2008I) . 



<wo 



200 



km s 1 (|Cerruti-Sola fc Perinottolll985t iGuerrero et al 



2010). In this interacting stellar- wind model (|Kwok 



19831 ). the physical structure of a PN is similar to that 
of a wind-blown bubble and will comprise (1) a cen- 
tral cavity filled with shocked fast wind at temperatures 
of 10 7 -10 8 K, (2) a dense shell of swept-up AGB wind 
at 10 4 K (the bright ring seen in optical images), and 
(3) an outer envelope of ionized AGB wind material re- 
shaped by the leading shock s e t up by ionization (cf. 
Schmidt- Voigt fc Hopped Il987t iMarten fc Schonberner 



19911 : lMellemalll994t iVillaver et al. 1 12003 iPerinotto et al 
2004). At the interface between the shocked fast wind 



and the swept-up AGB wind, a contact disco ntinuity 
forms and heat conduction is expected to occur (jSpitzerl 
1962). The resulting mass evaporation from the dense 
nebular shell into the hot interior lowers the tempera- 



nieves@iaa.es 

1 Institute- de Astroffsica de Andaluci'a, IAA-CSIC, c/ Glori- 
eta de la Astronorma s/n, 18008 Granada, Spain 

2 Department ol Astronomy, University of Illinois, 1002 West 
Green Street, Urbana, IL 61801, USA 

3 Leibniz-Institut fur Astrophysik Potsdam (AIP), An der 
Sternwarte 16, 14482 Potsdam, Germany 



Chandra and XMM-Newton observations of diffuse 
X-ray emission from PNe have been used to investi- 
gate the physical properties of the hot interior gas, 
which shows plasma temperatures of (1-3) xlO 6 K 
and X-ray luminosities Lx=2x 10 30 -3x 10 32 erg s" 1 
Ce.g. JChu et alJl2004t [Guerrero et al.ll2005HMontez et al.1 
[20051 : IKastner et al.l 120081 IRuiz et al.l 120111) ! The X- 
ray morphology, low plasma temperatures, and mod- 
erate Lx are quantitatively consistent with those ex- 
pected from bubble models with he at conduction 
(|Steffen. Schonberner. fc Warmuth I [20081 ) . 

To make rigorous comparisons between observations 
and model predictions of hot gas in PN interiors and 
to achieve statistical significance in the comparisons, the 
small sample of PNe with known diffuse X-ray emission 
needs to be enlarged. While it has been observed that 
PNe possessing detectable diffuse X-ray emission exhibit 
sharp inner shells of swept-up AGB wind, the reverse 
may not be true (jKastner et al.l[2~012t ). A more reliable 
diagnostic to ensure the detection of diffuse X-ray emis- 
sion from PNe is needed. To this end, we have utilized 
the collisionally ionized O VI expected in interface layers 
to diagnose the presence of hot gas. Using archival Far 
Ultraviolet Spectroscopic Explorer (FUSE; IMoos et al.l 
l2000f ) observations of PNe, we find three cases, IC418, 
NGC 2392, and NGC 6826, exhibiting nebular O vi ab- 
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TABLE 1 

Space Observations of IC418, NGC 2392, and NGC 6826 



FUSE UV Observations 


Object 


Program ID 


Aperture 


Date 


Processing Version 




'[sT 


IC 418 
NGC 2392 
NGC 6826 


LWRS 
LWRS 
MDRS 


P115U11 
B0320601 
P1930401 


2001 Dec 2 
2001 Feb 21 
2000 Aug 8 


CalFUSE v3.2.3 
CalFUSE v3.2.3 
CalFUSE v3.2.3 




4,440 
2,620 
5,800 






Chandra X-ray 


Observations 








Object 


Observation ID 


Instrument & Pointing 


Date 


Processing Version 


^-obs 

[ks] 


^cxp 

[ks] 


IC418 
NGC 2392 
NGC 6826 


7440 
7421 
7439 
8559 


ACIS-S S3 
ACIS-S S3 
ACIS-S S3 
ACIS-S S3 


2006 Dec 12 

2007 Nov 13 
2007 Jun 11 
2007 Jul 24 


DS 7.6.9 
DS 7.6.11 
DS 7.6.10 
DS 7.6.10 


51.1 
58.1 
34.5 
15.0 


50.4 
57.4 
34.1 
14.9 


HST Optical Observations 


Object 


Program ID 


Instrument & Pointing 


Date 


Filter 




[sf 


IC 418 
NGC 2392 
NGC 6826 


8773 
8499 
6117 


WFPC2-PC 

WFPC2-WF3 

WFPC2-PC 


2001 Oct 3 
2000 Jan 11 
1996 Jan 27 


F502N 
F656N 
F502N 




600 
100 
100 



sorption lines. We have subsequently obtained Chandra 
observations of these three PNe, and indeed diffuse X- 
ray emission is detected in all three cases. In this paper, 
we describe the FUSE observations of PNe in Section 
2, report the Chandra observations of the three PNe in 
Section 3, discuss the X-ray results in the framework of 
models of PNe formation and evolution in Section 4, and 
summarize the main conclusions in Section 5. 

2. O VI DIAGNOSTIC OF HOT GAS 

A PN with hot interior gas will possess an interface 
layer in which the temperature falls from 10 6 K near the 
hot gas to 10 4 K in the cool nebular shell. Thermal colli- 
sions in the interface layer produce highly ionized species, 
such as C IV, N v, and O VI, whose excitation poten- 
tials are 47.8, 77.5, and 113.9 eV, and whose fractional 
abundance s peak at — lxlO 5 , 2xl0 5 , a nd 3xl0 5 K, re- 
specti velv dShull fc van Steen bcrg 1982). Cloudy simula- 
tions ([Ferland et al.lll998T ) show, however, that such ions 
can also be produced by photoionization in appreciable 
amounts for luminous central stars of planetary nebu- 
lae (CSPNe) hotter than -35,000, 85,000, and 140,000 
K, respectively. Being the hardest to be produced by 
photoionization, O vi is thus the best choice among the 
three species for diagnosing interface layers and hot gas 
in PNe. O vi has a resonance doublet at AA1031. 9, 1037.6 
A, a primary spectral feature targeted by FUSE. 

The O vi doublet from the interface layer can be de- 
tected as emission features in FUSE spectra if the en- 
trance aperture does not include the CSPN, or as absorp- 
tion features if the entrance aperture includes the CSPN. 
FUSE observations of NGC 6543 and NGC 7009, two 
PNe with diffuse X-ray emission, indeed detected nebu- 
lar O vi emission from t he interface layers ( Grucnd l~et al.l 
I500i lining et^[2002h . 

The FUSE arc hive contains high-d ispersion spectra 
(A/AA w20,000, iSahnow et all l2000h of a large num- 
ber of CSPNe. These observations provide two pieces of 
information pertinent to the hot interior gas. First, the 
stellar O VI lines may show P Cygni profiles indicating 
the existence of a fast stellar wind that is needed to pro- 



duce the hot interior gas in a PN. Second, there may be 
narrow nebular O vi absorption features superposed on 
the stellar P Cygni profile, and if the CSPN is cooler than 
140,000 K, the nebular absorption must originate from 
an interface layer, requiring the existence of hot interior 
gas. FUSE makes simultaneous observations through the 
30"x30" low (LWRS), 4"x20" medium (MDRS), and 
l'.'25x20" high (HIRS) resolution apertures, that are off- 
set from one another. It is possible that some FUSE 
observations of CSPNe had one or two off-source aper- 
tures falling within the nebula, allowing us to search for 
O vi emission from the interface layer. 

We have used FUSE observations of the nebular shells 
of a dozen PNe to search for nebular O VI emission, and 
analyzed FUSE observations of ~60 CSPNe to search 
for narrow nebular O vi absorption lines in their spec- 
tra. Each observation is typically composed of multiple 
spectra. The individua l exposures were processed us- 
ing the CalFUSE v3.2.3 (jDixon et al.ll2007D to obtain an 
extracted spectrum. Prior to combining the individual 
spectra, the relative alignment of each was checked by 
performing a cross correlation. Typical offsets were less 
than a few km s _1 for the LiFIA module. The spec- 
tra were then combined using a weighted average after 
accounting for the small offsets in wavelength. 

Among the PNe with FUSE observations, IC418, 
NGC 2392, and NGC 6826 display narrow nebular O vi 
absorption lines superposed on the spectra of their 
CSPNe (Figure CD see Table Q] for the details of these 
observations). We have thus obtained Chandra X-ray 
observations of these three PNe. Note that the Chandra 
observation of NGC 2392 was a follow-up of an earlier 
XMM-Newton observation ([Guerrero et alJl2005f ). 

3. CHANDRA OBSERVATIONS 

The Chandra X-ray Observatory {CXO) observed 
IC418, NGC 2392, and NGC 6826 using the array for 
spectroscopy of the Advanced CCD Imaging Spectrom- 
eter (ACIS-S). All three PNe were imaged on the back- 
illuminated CCD S3 using the VFAINT mode. Details 
of the observations are given in Table [TJ The Chandra 
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Fig. 1. — FUSE spectra showing the region around the O VI lines for the central stars of IC 418 (top panels), NGC 2392 (middle panels), 
and NGC 6826 (bottom panels). For each object the upper panel shows a wide spectral range with the broad Ly/3 and O VI P Cygni profiles 
highlighted (darker grey). The shaded region in the upper panel is expanded to better show narrow absorption features superposed on the 
broader spectrum. Absorption from the H2 Lyman and Werner bands are marked in red and ionic transition are marked with blue. The 
narrow, nebular O VI components are also marked with vertical dashed lines. 



observations of NGC 6826 were split into two segments 



due to scheduling issues. 
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Fig. 2. — Chandra ACIS-S raw (left) and smoothed (center) X-ray images and HST optical narrow-band images (right) of IC418 (top 
row), NGC 2392 (middle row), and NGC 6826 (bottom row). The X-ray contours extracted from the smoothed X-ray images have been 
overlaid onto the optical images. The black and grey (only for IC418) contours mark the X-ray levels corresponding to 10%, 25%, 50%, 
75%, and 95% of the peak intensity of the diffuse emission. Two additional white contours defining the 25% and 75% of the peak intensity 
of the central source are shown for NGC 2392 and NGC 6826. 
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Fig. 3. — Chandra and HST composite pictures of IC418 (up- 
per), NGC2392 (center), and NGC 6826 (bottom). X-ray emission 
is shown in blue, and optical Ha and [N II] line emission in green 
and red, respectively. The field of view (FoV) of the pictures is 
Ri27" for IC 418, «63" for NGC 2392, and ^41" for NGC 6826. 



E (keV) 

Fig. 4. — Chandra ACIS-S S3 background-subtracted spectra 
of the diffuse X-ray emission of IC418 (upper), NGC 2392 (cen- 
ter), and NGC 6826 (bottom). The histograms correspond to mod- 
els that describe well the observed X-ray spectra of IC 418 and 
NGC 6826, while it is the best-fit model for NGC 2392. 



The background emission level was mostly stable for 
the observations of the three PNe, with NGC 2392 having 
a quiescent background level, and IC 418 and NGC 6826 
having abnormally high background levels. Besides a 
short spike of high background that affected the obser- 



vation of IC 418, no further period has been excised from 
the original data, and therefore the differences between 
the observation times i b s and the useful exposure times 
i exp in Table Q] are mostly associated with the removal 
of dead-time periods. All subsequent analyses have been 



6 



Ruiz et al. 



performed using the Chandra Interactive Analysis of Ob- 
servations (CIAO) software package version 4.1.2 and 
HEASARC XSPEC vl2.3.0 routines (|Arnaudlll996l) . 

The Chandra ACIS-S X-ray images of IC418, 
NGC 2392, and NGC 6826 are displayed in Figures [Hand 
[3J The raw X-ray images in Figure[5]have been extracted 
using the natural ACIS-S pixel size of 0'.'5. The smoothed 
images in Figures [2] and [3] have been processed using the 
CIAO task "csmooth" with a circular Gaussian kernel 
with size up to 4 pixels («2'.'0) and a fast-Fourier trans- 
form (FFT) convolution method. These images reveal 
the presence of soft diffuse X-ray emission in the three 
PNe. A comparison with archival HST WFPC2 narrow- 
band optical images selected to emphasize the innermost 
shells and small-scale nebular features (Tab. [T]) shows 
that this diffuse X-ray emission is confined within the 
innermost nebular shell (Figures [5] and [3]). We describe 
below in more detail the spatial and spectral properties 
of this diffuse emission. The analysis of X-ray point- 
sources at the CSPNe of NGC 2392 and NGC 6826 will 
be reported in an upcoming paper (Guerrero et al., in 
preparation) . 

3.1. IC'418 

The Chandra observations of IC 418 detect diffuse 
emission at an ACIS-S S3 count rate of 0.62±0.12 cnt 
ks" 1 in the energy range 0.3-2.0 keV for a total of 32±6 
counts. The X-ray emission is confined within the high 
excitation, innermost 3'.'5x5'/0 shell traced by the emis- 
sion in the [O in] line (Figs. [Hand [3]). Owing to the 
small angular size of the X-ray-emitting region and small 
number of counts, the spatial distribution of the X-ray 
emission cannot be analyzed in detail. Both the raw 
image and contours on the top panels of Figure [5] sug- 
gest that the diffuse X-ray emission is elongated along 
the major axis of the inner shell. The peak of the X-ray 
emission is located ~(X'5 towards the northwest of the lo- 
cation of the central star. The spatial distribution of the 
diffuse X-ray emission of IC418 relative to the optical 
shell and CSPN is r eminiscent of that of BD+30°3639 
(jKastner et al.ll2000h . 

The ACIS-S S3 spectrum of IC418 peaks at -0.7 keV 
(Figure [4]). A spectral fit is obviously not possible, but 
the observed spectrum can be reasonably well described 
by a thin plasma emission model with a plasma tem- 
perature of 0.26 keV and nebular chemical abundances 
absorbed by a colum n density of Njj ~ lxlO 21 cm" 2 
(|Pottasch et al.1120041 ). In this model, the observed X-ray 
flux in the 0.3-2.0 keV band is 2.5xl0" 15 erg cm" 2 s" 1 , 
and the intrinsic X-ray l uminosity i s 8.4xl0 29 erg s _1 
for a distance of 1.2 kpc (|Frewl l2008). These values are 
compiled in Table [2] which includes other X-ray proper- 
ties such as the hot bubble radius and unabsorbed X-ray 
surface brightness. 

3.2. NGC 2392 

The Chandra observations of NGC 2392 resolve the X- 
ray emis sion detected by pre vious XMM-Newton obser- 
vations (jGuerrero et al.ll2"005[ ) into a point-source at the 
central star and diffuse emission within the innermost 
shell of this nebula (Figures [5] and [3J . The diffuse emis- 
sion is well confined within the 15'.'2x 17" 6 innermost 
nebular shell of NGC 2392. Its ACIS-S S3 count rate 



in the energy range 0.3-1.5 keV is 9.3±0.4 cnt ks -1 , and 
a total of 530±25 counts are detected. The diffuse X-ray 
emission does not show a limb-brightened morphology, 
but it is brighter in the central region and towards a 
northern region inside the central cavity. The spatial 
distribution of this X-ray emission is roughly consistent 
with that expected for an ellipsoidal shell filled with hot 
gas. 

The ACIS-S S3 spectrum of the diffuse emission of 
NGC 2392 peaks at 0.5-0.6 keV and shows a plateau 
of fainter emission between 0.7 and 1.0 keV. Very lit- 
tle emission is detected above 1.0 keV. The ACIS-S 
spectrum can be well fit by a thin plasma MEKAL 
emissi on model. The an alysis of the XMM-Newton 
data (|Guerrero et al.l I2005D suggested that the X-ray- 
emitting plasma in NGC 2392 had N/O and Ne/O abun- 
dance ratios gre ater than the respective nebular values 
of 0. 4 and 0.2 (jBarkerl fl99lt IHenrv. K witter, fc Bates! 
2000). The more recent nebular abundance study of 
NGC 2392 bv lPottascheTaLl (120081) suggests N/O=0.65 
and Ne/O=0.30. These nebular abundances, together 
with the absorption column density Ah=9x10 20 cm" 2 
converted from the optical extinction, provide a reason- 
ably good fit to the observed X-ray spectrum (reduced 
X 2 =1.12) for a plasma temperature of fcT=0.18±0.04 
keV. In this model, the observed flux in the 0.2-1.5 keV 
band is 3.9xl0" 14 erg cm" 2 s _1 , and the intrinsic X-ray 
luminosity i s 1.8 xlO 31 erg s _1 for a distance of 1.28 kpc 
(|Prewl [20081 ) . Other X-ray properties of NGC 2392 are 
listed in Tabled 

The comparison between the current spectral fit and 
that performed on the XMM-Newton data obtained 
on 2004 April 2 needs to take into account that 
XMM-Newton did not resolve the point source at the 
central star from the diffuse emission. The physi- 
cal conditions implied by the spectral fits to the dif- 
fuse emission detected by Chandra (A^h=9x10 20 cm" 2 , 
fcT=0.18±0.04 keV) and by XMM-Newton (7V H =8xl0 20 

cm" 2 , fcT=0.18t°;o3 keV ) are consistent, but the 
observed diffuse X-ray flux derived from Chandra 
(/ x =3.9xl0" 14 erg cm" 2 s" 1 ) is -35% lower than that 
from the XMM-Newton (/x=6.0x 10" 14 erg cm" 2 s" 1 ). 
Part of this difference is caused by the contribution of 
the central point source, of which the observed flux is 
1.4x10" 14 erg cm" 2 s" 1 (Guerrero et al., in prepara- 
tion). The remaining ~12% difference is only slightly 
larger than the calibration uncertainties between Chan- 
dra ACIS and XMM-Newton EPIC0. 

3.3. NGC 6826 

The Chandra observations of NGC 6826 detect a point- 
source at its central star (Guerrero et al., in preparation) 
and diffuse emission that, as is the case with IC418 and 
NGC 2392, is confined within the ~12'.'4x7'/4 innermost 
shell of the nebula (Figs. [5] and [3]). The diffuse emis- 
sion is detected at an ACIS-S S3 count rate of 1.97±0.21 
cnt ks -1 in the energy range 0.3-2.0 keV for a total of 
96±10 counts. The spatial distribution of the diffuse X- 
ray emission shown by the contours on the bottom panels 

4 As de scribed by the XMM-Newton Calibration Team in 
|http: / /xmm.esac.esa.in t / external / xmm_calibration| these differ- 
ences can amount to 1U%, especially m the sott energy range below 
1 keV. 
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TABLE 2 

X-ray Detections of Diffuse Emission from PNe with Bubble Morphology 



Object 


Telescope & 
Instrument 


d a 


Hot Bubble 
Radius 


Count Rate 


fx 




Sx b 




Lx h 


kTx 


Ref. 






[kpc] 


[PC] 


[cnt ks" 1 ] 


[erg cm - 2 




[ergcm 2 s ^^arcsec 2 ] 


[ergs" 1 ] 


[keV] 




BD+30°3639 


CXO ACIS-S 


1.30 


0.013 


244±4 


5.7x10" 


■13 


2.2x10" 


■14 


2.7xl0 32 


0.233±0.007 


1 


IC418 


CXO ACIS-S 


1.20 


0.012 


0.62±0.12 


2.5x10" 


-15 


8.4x10" 


-17 


8.4xl0 29 


0.26 


2 


NGC40 


CXO ACIS-S 


1.02 


0.10 


2.8±0.9 


1.3x10" 


-14 


8.3x10" 


-17 


2.1xl0 31 


0.09 


3,4 C 


NGC 2392 


CXO ACIS-S 


1.28 


0.050 


9.3±0.4 


3.9x10' 


-14 


1.0x10" 


-16 


1.8xl0 31 


0.18±0.04 


2 


NGC 3242 


XMM EPIC 


1.0 


0.041 


31.3±1.6 


4.2x10" 


-14 


5.8x10" 


-17 


7.3xl0 30 


0.202±0.011 


5 


NGC 5315 


CXO ACIS-S 


2.62 


0.017 


12.4±0.7 


1.0x10" 


-13 


1.5x10" 


-14 


2.9xl0 32 


0.224±0.022 


■1 


NGC 6543 


CXO ACIS-S 


1.50 


0.037 


42.4±0.9 


9.5x10" 


14 


5.7x10" 


-16 


6.5xl0 31 


0.145±0.010 


6 


NGC 6826 


CXO ACIS-S 


1.30 


0.028 


1.97±0.21 


9.0x10" 


■15 


3.4x10" 


-17 


2.0xl0 30 


0.20 


2 


NGC 7009 


XMM EPIC 


1.45 


0.051 


61.5±1.7 


7.2x10" 


14 


1.9x10" 


-16 


4.4xl0 31 


0.152±0.015 


7 


NGC 7027 


CXO ACIS-S 


0.89 


0.012 


14.0±0.9 


3.1x10" 


-14 


5.2x10" 


-14 


1.3xl0 32 


0.26 
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Ref eren ces. — (1) IKastner et al.1 ifeoOOl). (2) th i s wor k, (3) IMontez et all 120051) . 14) IKastner et all 120081) . (5) IRuiz et al.l 120111) . 16) IChu et al.1 
12001fl . (7) IGuerrero et all 120051) . (8) IKastner et all 120011) . 

a Distances adopted from IMellemal 12004D and IFrewi 120081) . 

k Unabsorbed X-ray surface brightness and luminosity in the energy range 0.3-2.0 keV. 

c Montcz ct al. (2005) mistakenly reported the absorbed X-ray luminosity of NGC 40 whose value was later corrected by IKastner et al.1 | [20ul) . 



of Figure [2] hint at a patchy distribution which seems to 
be better described by a centrally filled shell than by the 
limb-brightened morphology that would produce a thin 
shell. 

The ACIS-S S3 spectrum of NGC 6826 is soft, with 
a plateau of emission between 0.35 and 0.7 keV over 
which a subtle peak at ^0.45 keV can be seen. A 
weaker emission peak is detected at ~0.9 keV, and no 
noticeable emission is detected above 1.0 keV. The spec- 
tral shape is suggestive of emission from an optically 
thin plasma, but, as with IC418, no reliable spectral 
fit is possible due to the small number of counts. As- 
suming the chemical abundances and absorbing column 
density of Ah ~ 1 x 10 20 cm" 2 derived for NGC 6826 
(jSurendiranath fc Pottascr] 120081 ). the spectrum of its 
diffuse X-rays can be reasonably well described by a 
plasma emission model for a plasma temperature of kT 
= 0.2 keV and the observed X-ray flux is 9.0xl0~ 15 
erg cm -2 s _1 in the 0.3-2.0 keV energy band. The in- 
trinsic X-ray luminosity for this model and energy ban d 
is 2.0xl0 30 erg s" 1 for a distance of 1.3 kpc (|Frewil2008[) . 
Further X-ray properties of NGC 6826 are listed in Ta- 
bleU 

4. DISCUSSION 

The background-subtracted spectra of IC418, 
NGC 2392 and NGC 6826 have been described using 
optically-thin thermal plasma models with nebular 
abundances and their basic parameters (temperature 
and X-ray flux and luminosity in the energy band 0.3-2.0 
keV) determined in previous sections. The detection of 
diffuse X-ray emission from these three PNe, together 
with NGC 6543 and NGC 7009 (IGruendl et all T2001 
llping et ail |2002[ ). testifies that nebular O VI emis- 
sion and/or absorption is an excellent diagnostic for 
the presence of hot bubbles in PNe with sharp shell 
morphology. 

Similarly, the detection of diffuse X-ray emission in 
a PN can be used to forecast the presence of nar- 
row O VI absorptions in the spectrum of its CSPN. 
Am ong the PNe with d iffuse X-ray emission (this pa- 
per; IKastner et all [20121 ) . there are another four with 



available FUSE observations of their central stars 
in the spectral range of the O VI lines, namely 
BD+30°3639, NGC 40, NG C 2371-2, and NGC 7662 
((Guerrero fe; De Marcol2 013'). A careful scrutiny of these 
spectra reveals narrow O VI absorptions of the A1031.9 A 
line in NGC 40, NGC 2371-2, and possibly in NGC 7662, 
as well as a possible absorption of the A1037.6 A line in 
NGC 2371-2. As for BD+30°3639, the search for narrow 
O VI absorptions in its FUSE spectrum is inconclusive 
because the large number of atomic and H2 absorption 
lines that dominate its stellar continuum at this wave- 
length range. To sum up, 8 out of the 9 PNe with diffuse 
X-ray emission also show O vi narrow absorptions in the 
stellar continuum of their CSPNe. 

The firm correlation between diffuse X-ray emission 
and O Vi narrow absorptions in the stellar continuum 
of PNe with sharp shell morphology provides strong evi- 
dence for a conduction layer between the hot interior and 
the cool nebular shell of PNe. The physical structure 
(how the density and temperature vary with radial dis- 
tance) of this conduction layer and the amount of highly 
ionized species present at this interface (0 +5 , N +4 ) de- 
pend on the efficiency of thermal conduction, although 
the O VI luminosity seems rather insensitive to those ef- 
fects (|Steffen. Schonberner. fc Warmuth 1 120081 ). We ex- 
pect our awarded HST STIS observations of the N v 
and C iv line emission from the interfaces in NGC 6543 
(PI: M.A. Guerrero) and in the Wolf-Rayet wind-blown 
bubble S 308 (PI: Y.-H. Chu) would help us obtain more 
information on the physical structure of conduction lay- 
ers. 

4.1. Hot Bubbles of Planetary Nebulae 

Chandra and XMM-Newton have yielded a number of 
detections of diffuse X-ray emission confined within the 
innermost closed shells of PNe, the so-called hot bubbles. 
Sources with available X-ray luminosities and tempera- 
tures in the liter atur^ are summarized in Table [5] where 
we have excluded the bipolar PNe Mz 3 and NGC 7026 

5 Recently, Kastncr ct al. (2012) has reported the detection 
of diffuse X-ray emission in the elliptical PNe NGC 2371-2 and 
NGC 7662, but no estimates of their X-ray temperatures and lu- 
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TABLE 3 

Properties of the stellar wind and diffuse X-ray emission of PNe 



Object 


Spectral 


T e fi tioo 


logL/L a M a 


log(L x /L.) b log(L x /L w ind) b T shock c T x /T shock Ref. 




Type 












[kK] [km s- 1 ] 


[Mo yr- 1 ] 


[10 6 K] 



BD+30°3639 


[WC9] 


47 


700 


3.77 


1.8xl0" 6 


-4.93 


-3.01 


14.8 


0.2 


1 


IC418 


Of(H) 


39 


700 


3.72 


3.3xl0~ 8 


-7.39 


-3.78 


7.4 


0.4 


2 


NGC40 


[WC8] 


71 


1000 


3.42 


l.lxlO" 6 


-5.69 


-4.22 


30.2 


0.035 


1 


NGC 2392 


06f 


45 


300 


3.82 


8.4xl0~ 9 


-6.16 


-1.12 


1.4 


1.5 


3 


NGC 3242 


O(H) 


75 


2400 


3.42 


3.4xl0~ 9 


-6.15 


-2.93 


87.0 


0.03 


2 


NGC 5315 


[W04] 


76 


2400 


3.74 


1.6xl0~ 6 


-4.87 


-4.00 


174.0 


0.015 


1 


NGC 6543 


wels 


60 


1450 


3.54 


3.3xl0- g 


-5.32 


-2.53 


31.8 


0.05 


2,3 


NGC 6826 


03f(H) 


14 


1200 


3.81 


4.7xl0" 8 


-7.10 


-4.03 


22.0 


0.1 


2 


NGC 7009 


0(H) 


87 


2770 


3.56 


1.2xl0" 9 


-5.51 


-1.82 


115.9 


0.015 


4,5,6 


NGC 7027 


198 




3.89 




-5.37 
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References. (1) [Marcolino ct al. 120071), (2) IPauldrach, et al.l J200l) , (3) IHerald & Blanch] 120111) , (4) IMendez et~aT1 11988ft , (5) 

ICerruti-Sola fe Perinottd Jl989h . f6) |Iping etal] j2002fl , (7) |Latter et alJ j2000ft . 

a Stellar luminosities and mass- loss rates adopted from the references in the last column and scaled with the distances given in Table [2] 
k X-ray luminosity in the 0.3-2.0 keV energy band. 

c The shock temperature of the [WR] stars BD+30°3639, NGC 40, and NGC 5315 is at given wind velocity about twice that of normal CSPNe, 
according to the larger mean molecular weight of their stellar winds. 



(IGuerrero. Chu. fe Mirandal l200i IGruendl et all [20061 
iClark et al.l 12012ft . Table l2l compiles the X -ray proper- 
ties of the PNe with bubble morphology, including their 
observed X-ray flux (fx) an d intrinsic surface bright- 
ness (Sx) and luminosity in the ener gy band 0.3- 2.0 keV, 
scaled to the distance determined bv lFrewl ()2008ft . as well 
as the X-ray temperature and hot bubble radiusQ. 

It is interesting to compare the physical properties of 
these CSPNe and their winds with the physical parame- 
ters of hot gas inferred from the diffuse X-ray emission. 
The properties of stellar winds of these PNe have been 
compiled in Table[3]that includes the stellar spectral type 
(column 2), effective temperature (column 3), wind ter- 
minal velocity (column 4), stellar luminosity (column 5), 
and mass-loss rate (column 6) as provided by the refer- 
ences listed in the last column of the table. An inspection 
of the spectral type of the CSPNe reveals that this sam- 
ple consists of three H-poor [WR] CSPNe (BD+30°3639, 
NGC 40, and NGC 5315), and six H-rich CSPNe (IC418, 
NGC 2392, NGC 3242, NGC 6543, NGC 6826, NGC 7009, 
and NGC 7027). The wind terminal velocities, stellar lu- 
minosities, and mass-loss rates have been derived from 
model atmosphere analyse s of UV and opti cal stellar 
lines, except for NGC 7027 (|Latter et al.l l2000). We note 
that the stellar luminosities and mass-loss rates provided 
here have been scaled with respect to those given in the 
original references by the distance given in the third col- 
umn of Tabled 

The X-ray luminosities and temperatures of these neb- 
ulae are then compared to their stellar (L*) and wind 
mechanical (£ W md = \ M v^) luminosities in columns 
7 and 8 of Table |3] Finally, the post-shock temperature 
of the stellar wind (T s h oc k), estimated assuming an adia- 



minosities are available. Similarly, diffuse X-ray emission has been 
detected within the innermost shell of IC4593 (Guerrero et al., in 
preparation). 

6 We note here that the values o f the hot bubb l e radi i of the 
PNe in Table[2]that were analyzed by Kastncr ct al. (2008) clearly 
differ with their values because they used diameters but quoted 
them as radii. 



batic shock as 

rp _ Ml! _ 

J- shock — , V,, 

16 k 



ee 2.3 x 10V 



1000 km s~ 



[K], 



where fj, is the mean molecular weiglrfl, is listed in col- 
umn 9 and compared with the plasma temperature de- 
rived from X-ray spectra in column 10. 

These observationally determined properties of hot 
plasma in PNe are compared in Figure [5] to theoreti- 
cal predictions based on ID hydro dynamical simulations 
performed with the code NEBEL (fPerinotto et al.l 1200 
Steffen. Schonberner. k Warmuth I 120081 : iSteffen et al 
2012ft . These simulations take into account the en 



ergy transfer due to heat conduction and the evolu- 
tion of stellar mass-loss rate and terminal wind veloc- 
ity that determine the mechanical energy input into 
the hot bubble. The plots presented in Figure [5] 
are thus an increment and refinement o f the plots in 
ISteffen. Schonberner. k Warmuth I (|2008ft for a larger 
sample of PNe with consistently revised stellar properties 
and diffuse X-ray parameters. 

Before starting a discussion on the location of the PNe 
in our sample on these plots, we would like to note that 
the models presented in Figure [S] have been developed 
exclusively for PNe with H-rich central stars. Therefore, 
the PNe with [WR] central stars BD+30°3639, NGC 40, 
and NGC 5315 (marked in these plots by "pink" symbols) 
are not expected to follow these tracks. Indeed, the data 
points of these sources appear notably far from these the- 
oretical tracks, e.g., BD+30°3639 and NGC 5315, with 
young [WR] CSPNe, exhibit the highest values of Lx 
(Fig. [5la) and largest Lx/L± ratios (Fig. [5lb, d, and 
e). Interestingly, the X-ray temperatures predicted by 
the models for H-rich stellar winds do not differ criti- 
cally from the temperature of the hot plasma derived 
for PNe with [WR] CSPNe (Fig. Of), despite that they 
have more massive stellar winds with higher mechanical 

7 The mean molecular weight for fully ionized H-rich stellar 
winds can be approximated as Ri0.6, whereas for fully ionized H- 
poor stellar winds is Ril.2. 
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Fig. 5. — X-ray luminosity (Lx) of the hot bubble of PNe for the energy range 0.3-2.0 keV (6.2-41 A) and X-ray temperature (Tx) 
as functions of T c fj, the stellar effective temperature (panels a, c, and e), X-ray temperature (panel e), and V w i n d, the wind terminal 
velocity (panels d and f). In panels b, c,d, and e, Lx has been scaled with the stellar luminosity (L*) and wind mechanical luminosity 
(L wind ) to remove the dependence on distance. The tracks correspond to the theoretical models HC2 with heat conduction for CSPNe 
with masses 0.565 Afp (green), . 595 M fl (red), and 0.696 Mq (blue) according to the prescriptions of the second method used by 
Steffcn, Schonbcrncr, & Warmuth (2008), i.e., the heat flux does not exceed the saturation limit. The dashed line in panel f represents the 
post-shock temperature of an adiabatic shock for the case of H-rich stars, i.e., T s h oc k described by Eq. 1 for /i=0.6. Black dots correspond 
to objects with a H-rich CSPN, whereas pink dots refer to sources with a [WR]-type CSPN. The error-bars in the data points correspond 
to the 1-<T uncertainty in our determination of Lx and Tx (the uncertainty in distance is not included in the error- bar of Lx)- Typical 
uncertainties of 0.15 dex for L* and 0.8 dex for L w ; n( j have been convolved with the error-bars of Lx for the Lx/L* and Lx/L w ; n( j ratios. 
Similarly, typical uncertainties of 0.05 dex in T c ff and 100 km s _1 in V„i n( j are presented. 



luminosities than their H-rich counterparts. Theoretical 
work is underway to investigate t he X-ray properties of 
hot bubbles around [WR] CSPNe (ISteffen et al. l lMl) . 

Overall, there is a good agreement between the global 
properties of the hot gas in PNe with H-rich stars 
("black" symbols in Figure [SJ and the model predic- 
tions, i.e., most of the data points are bracketed within 
the tracks predicted by the models. The most notable 



differences appear in the comparison between models 
and observations when the stellar wind mechanical lu- 
minosity, L W ind, is considered (Fig. [5]-c). The disper- 
sion of the data points in this figure can stem from 
the large uncertainties in the mass-loss estimates, but 
also from the comparison of the present-time mass loss 
with X-ray properties that depend on the mass-loss rate 
integrated over time. The difficulties for such com- 
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Fig. 6. — Comparisons of the X-ray plasma temperature with hot 
bubble radius of PNe (top) and CSPN effective temperature (mid- 
dle), and X-ray surface brightness with CSPN effective temperature 
(bottom). The theoretical tracks of the variations of Tx and Sx im- 
ply that the observations comply with the expectations of the theo- 
retical models HC2 including heat conduction for central stars with 
masses 0.565 Mq (green), 0.59 5 Mq (red), 0.625 Mq (orange), and 
0.696 Mq (blue) adopted from Stcffcn. Schonberncr. & Warmuth 
(2008). Error- bars as in Figure [5] 



parison are well illustrated b y NGC 7009, whose low 
mass- loss rate, as provided by ICerruti- Sola fc Perinottol 
1989) , is increased by o ther aut hors by factors of ~1. 8 
Tinkler fc L amersl 12001 . ~3.6 dBombeck et all 119861). 

On 



and ~10.7 (jCerruti-Sola fc Perinottol 119851 ). On the 
other hand, the increase of X-ray luminosity with time, 
as revealed by the evolution dependent stellar parame- 
ters T e ff and Voo shown in panels a, d, and e of Figure [SJ 
is consistent with the model predictions. Similarly, the 
little variation of Tx with the evolution predicted by the 
models is supported by the data points shown in panels 



b and f of Figure [S] 

As for the PNe with new data presented in this pa- 
per, we note that IC 418 and NGC 6826 seem to have 
low values of Lx and Lx/L+, and very notably small 
Lx/L W i n d (Fig. 03-c). We should emphasize that the X- 
ray luminosities of these two PNe have not been derived 
from spectral fits, but from a rough estimate of their 
X-ray temperatures and count rates due to the small 
number of detected counts (~32 cnts and ~96 cnts, re- 
spectively). The uncertainty in the X-ray temperature 
of these sources can amount up to 15-25%, whereas the 
count rate uncertainty is ~20% for IC418 and ~10% 
for NGC 6826. Taking into account the uncertainties in 
count rates, plasma temperatures, and absorption col- 
umn densities, the X-ray luminosities of these two sources 
shown in Table [2] can be higher by 50%. Even then, the 
locations of these two sources in the £x/£wind vs - ^eff 
(Fig. Ob) and Lx/L wind vs. T x (Fig. EJc) plots still fall 
below the theoretical tracks. Contrary to IC418 and 
NGC 6826, the data point of NGC 2392 is well above the 
theoretical tracks in the £x/£wind vs. Tx plot (Fig.Oc). 
We will discuss later in section 4.2 the X-ray peculiarities 
of NGC 2392. 

It has been suggested that the time-scales for the ef- 
fects of ste llar wind shock-heated plasma in PNe is short, 
<5,000 yr (|Ruiz et al.ll20TTI: iKastner et al.|[2012l ). To in- 
vestigate the time evolution of the hot bubbles of PNe, we 
have plot in Figure [6] the X-ray temperature of the PNe 
in our sample with H-rich stellar winds against the radius 
of their hot bubbles (top panel) and the effective temper- 
ature of their CSPNe (middle panel) . The top panel of 
this figure shows an apparent anti-correlation between 
the X-ray temperature and hot bubble radius, i.e., as 
the hot bubble expands, the X-ray temperature of the 
plasma seems to decrease. The comparison of the data 
points with the theoretical tracks on this plot reveals 
that we may be comparing PNe descending from progen- 
itors of different initial masses. This is more clearly illus- 
trated in the middle panel plotting Tx against T e g; PNe 
with more massive central stars produce hotter X-ray- 
emitting plasma because they reach a high wind power 
in a short time-scale, when the hot bubble is still small in 
size. The plot of X-ray surface brightness (Sx) against 
the star effective temperature, in the bottom panel of 
Figure [6l further supports this conclusion. The value of 
Sx does not vary much as the star moves along the hor- 
izontal track of the post-AGB evolution, but it depends 
strongly on the mass of the central star. PNe with mas- 
sive central stars rapidly inject larger amounts of wind 
mechanical energy into small-sized bubbles, resulting not 
only in higher plasma temperatures, but also in higher 
X-ray surface brightnesses. 

Finally, it is interesting to note the persistent absence 
of limb-brightening in the spatial distribution of dif- 
fuse X-ray emission from hot bubbles of PNe (see also 
IKastner et al.l 120121) . In particular, the diffuse X-ray 
emission from PC 418, NGC 2392, and NGC 6826 does 
not show clear evidence for limb-brightening (Fig. [2]), 
suggesting that the X-ray-emitting plasma may fill a 
significant fraction of the nebular shell. Whereas mod- 
els of hot bubbles with heat conduction predict a limb- 
brightened morphology for this diffuse emission, it has 
been noted that even a small amount of interstellar 
extinction may reduce significantly the center-to-limb 
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emission contrast because the soft emission from the 
cooler plasma close to the nebular rim is more easily ab- 
sorbed th an the harder emission from the hotter p lasma 
inwards (jSteffen. S chonbcr ner. fc Warmuthl l2008| ) . In 
this respect, the spatial distribution of diffuse X-ray 
emission from PNe differs notably from that observed 
in wind-blown bubbles around Wolf-Rayet stars, which 
show distinct lim b-brightened morphologies (e.g., S308, 
iToala et al.ll2012T) . 

4.2. NGC2392 - Over-luminous for Its Wind 

An inspection of the location of the PNe whose Chan- 
dra observations are presented in this paper on the dif- 
ferent plots shown in Fig. [5] suggests that IC 418 and 
NGC 6826 are generally a bit under-luminous, whereas 
NGC 2392 follows the theoretical tracks with the notable 
exception of its high Lx/^wind in Fig- EJ-c Such pecu- 
liar behavior of NGC 2392 is confirmed by comparing its 
values of £x/£ w ind an d ^x/^shock to those of the other 
PNe in Table [3] The diffuse X-ray luminosities of these 
PNe comprise just minute fractions of their stellar wind 
mechanical luminosities; however, the diffuse X-ray lu- 
minosity of NGC 2392 constitutes a significant fraction, 
^10%, of its wind mechanical luminosity. In absolute 
terms, the X-ray luminosity of NGC 2392 exceeds that 
of many PNe in Table [2] which is puzzling as the stellar 
wind of the CSPN of NGC 2392 has a relatively low mass- 
loss rate and a meager terminal velocity of 300 km s . 

The plasma temperature of NGC 2392 is also puzzling. 
While the plasma temperatures implied by X-ray spec- 
tra of all PNe in Table [2] are lower than the post-shock 
temperatures of their stellar winds, the observed plasma 
temperature of NGC 2392 is higher than the post-shock 
temperature expected from its 300 km s _1 wind terminal 
velocity. The influence of heat conduction is expected to 
lower significantly the plasma temperature from the post- 
shock temperature (dashed line in Fig. [5}f), which rises 
with increasing wind velocity, as it is certainly the case 
for NGC 6826. However, at low wind speeds, early in the 
evolution of the hot bubble, the influence of heat conduc- 
tion on plasma temperature is still small, and therefore 
the plasma temperature is not expected to differ substan- 
tially from the post-shock temperature of an adiabatic 
shock. Yet, the plasma temperature should not exceed 
the post-shock temperature. 

We conclude that the stellar winds of the PNe in our 
sample are able to power their X-ray-emitting hot gas 
except for NGC 2392. A wind terminal velocity higher 
than 300 km s^ 1 is needed for NGC 2392 to raise its ex- 
pected for Lx and Tx- We note that larger values of 
the terminal velocity of t his w ind have been r e porte d 
by IPauldrach et al.l (120041) a nd IKudritzki et~aTl (|2006l) . 
but Herald fc Bianchil (|20Tl have examined these val- 
ues and found them not compatible with the P-Cygni 
profiles seen in FUSE observations. Considering that an 
unusual 20 km s" 1 bipolar outflo w has been reported in 
NGC 2392 (jGieseking et al.l 119851 ) and that the nebular 
ionization and expansion velocity are anomalously high, 
it is possible that the CSPN has a hidden companion 
and their binary interactional contribute to the energet- 

8 Such interaction s exclude the faint companion reported by 
ICiardullo etaE! £1999) to be in a possible binary association be- 
cause at the distance of 2'.'65 from the CSPN of NGC 2392 it would 



ics of the nebular interi or (iDanehkar et al.ll2012l Mendez 
et al., in preparation). iGuerrero et a l. (2005) suggested 
that a fraction of the diffuse X-rays from NGC 2392 was 
associated with this bipolar outflow, an idea later pur- 
sued theoretically by lAkashi et al. (2008), but the spa- 
tial distribution of the diffuse X-ray emission revealed by 
Chandra does not support this hypothesis. Future inves- 
tigations of the CSPN of NGC 2392 are needed to solve 
the puzzle of the high X-ray luminosity of its hot interior. 

5. SUMMARY AND CONCLUSIONS 

We have obtained Chandra X-ray observations of three 
PNe, IC 418, NGC 2392, and NGC 6826, that display nar- 
row nebular O vi absorption lines superposed on the 
FUSE spectra of their CSPNe. Diffuse X-ray emis- 
sion is detected in each of these three PNe within the 
sharp innermost optical shell. These detections, to- 
gether with those in NGC 40, NGC 2371-2, NGC 6543, 
NGC 7009, and NGC 7662, amount to eight known PNe 
with hot (~ 10 6 K) interior gas in contact with the 
cool (~ 10 4 K) nebular shell whose conduction layers 
have been confirmed by the presence of collisionally pro- 
duced O vi. The physical structure of these PNe is 
thus consistent with that expected from bubble mod- 
els for H-rich stellar winds including he at conduction 
(Stcff en. Schonberner. fc Warmuth II2008D which predict 
a conduction layer with a steep temperature gradient be- 
tween the hot interior and the cool nebular shell. 

These models also have specific predictions as for the 
evolution in time of the global X-ray properties of bub- 
ble models. To test these predictions, we have compiled 
relevant information on the X-ray, stellar, and nebular 
properties of PNe with a bubble morphology. The ex- 
pectations of bubble models including heat conduction 
compare positively with the present X-ray observations 
of PNe with H-rich stellar winds, but those with H-poor 
[WR] stars present notable discrepancies. There is an ap- 
parent anti-correlation between X-ray temperature and 
hot bubble radius, but we note that this anti-correlation 
is cause by differences in stellar masses rather than evo- 
lutionary stages: PNe with massive central stars are ex- 
pected to produce hotter X-ray-emitting plasma inside 
small-sized hot bubbles of higher X-ray surface bright- 
ness. On the other hand, PNe with less massive central 
stars have a slower evolution during the post-AGB phase 
and inject less amounts of mechanical wind luminosities 
into larger hot bubbles, resulting in cooler X-ray-emitting 
plasma. 

Finally, we note that the stellar wind of PNe with bub- 
ble morphology and hot gas confined in their interiors 
seems capable of powering their X-ray emission, except 
for the notable case of NGC 2392. The low speed and 
small mass-loss rate of its stellar wind result in a low 
wind mechanical luminosity both in absolute and rela- 
tive terms compared to other PNe. We suggest that a 
binary companion may exist and contribute to high en- 
ergy processes that have provided the additional power 
to its hot gas. 
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